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Abstract 

We investigate the fractal structure of 2d quantum gravity, both for pure gravity 
and for gravity coupled to multiple gaussian fields and for gravity coupled to Ising 
spins. The roughness of the surfaces is described in terms of baby universes and 
using numerical simulations we measure their distribution which is related to the 
string susceptibility exponent ^string- 
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1 Introduction 



The fractal and selfsimilar structure of 2d quantum gravity is related to the entropy- 
(or string susceptibility-) exponent 7. This has been discussed in a recent paper |I[ 
where the structure of so-called baby universes was analyzed. It is convenient in the 
following discussion to consider 2d quantum gravity with an ultraviolet cut-off and 
we will consider the surfaces entering in the path integral as triangulated surfaces 
built out of equilateral triangles f||, |, U . In the case of surfaces of spherical topology 
a closed, non-intersecting loop along the links will separate the surface in two parts. 
The smallest such loop will be of length 3. It will split the surface in two parts. If the 
smallest part is different from a single triangle we call it (following the notation of 
|TJ) a "minimum neck baby universe", abbreviated "mimbu". The smallest possible 
area of a mimbu is 3 and the largest possible area will be Nt/2, where Nt is the 
number of triangles constituting the surface. 

In the case of pure 2d quantum gravity it is known that the number of distinct 
surfaces of genus zero made out of Nt triangles has the following asymptotic form 

Z(N T ) ~ e^ NT N}' 3 (1.1) 

where 7 = —1/2. For the models which can be solved explicitly and where c < 1 we 
have the following partition function: 

Z(ri = J2Z(N T )e-^ (1.2) 

Nt 

where Z(Nt) for large Nt is of the form ( |1.1| ), just with a different 7 = 7(c). 
For c = 1 it is known that there are logarithmical corrections to flLip , while the 
asymptotic form of Z(Nt) is unknown for c > 1, although it can be proven that it 
is exponentially bounded (@). If we assume ( |1.1| ) one can prove that the average 
number of mimbu's of area B on a closed surface of spherical topology and with 
area Nt (we use the notation area = ^triangles) is given by 

n NT (B) ~ (Nt - S) 7_2 S 7-2 (1.3) 

provided Nt and B are large enough. 

The above formalism is well suited for numerical simulations. The measurement of 
the exponent 7 has always been somewhat difficult. The first attempts used a grand 
canonical updating ([§, [7], |[), which generated directly the distribution (|1.2j). The 
disadvantage is that one has to fine-tune the value of /x to fi c . Later improved versions 
allowed one to avoid this , but one still had to perform independent Monte Carlo 
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simulations for a whole range of Nt and 7 still appeared as a subleading correction 
to the determination of the critical point \i c . These disadvantages disappear when 



keeps Nt fixed, and still in a single Monte Carlo simulation get a measurement 
of the distribution of mimbu's all the way up to Nt/2. One can therefore take 
a large Nt and make one very long run. This allows us to avoid the problems 
with long thermalization time. In addition the actual measurement of the mimbu 
distribution is easy. For a given thermalized configuration one has to identify all 
possible mimbu's associated with the configuration. This is done by picking up one 
link, Iq, and checking whether any links which have a vertex in common with Iq have 
a vertex in common which do not belong to Iq. This being the case we will have a 
minimal neck of length 3. For a given l there will always be two such, corresponding 
to the two triangles sharing Iq. But there might be additional ones and they will 
divide the surface into a mimbu and its "mother". By scanning over / 's, avoiding 
double counting and repeating the process for independent configurations we can 
construct the distribution of mimbu's. 

In the rest of this paper we report on the results of such numerical simulations. 

2 Numerical simulations 

2.1 Pure gravity 

The simulations were done on lattices of size ranging from 1000 to 4000 triangles 
Nt and we used the standard "link flip" algorithm [[| to update the geometry. We 
used of the order of 10 7 sweeps, where each sweep consists of Nt link flips. After 
thermalization we measured for each 10 th sweep the distribution of mimbu's, that 
is we counted all areas B > 1 enclosed by boundaries of length 3. The reason for 
performing the measurements so often is simply that they are not time consuming 
(the time it takes to make one measurement is comparable to the time it takes to 
perform one sweep). The distributions are shown in fig. 1. In order to extract 7 the 
distributions are fitted to equation (|1.3|). But as eq. ( J1.3| ) is only asymptotically 
correct deviations can be expected for small B. Thus a lower cut-off Bq has to be 
introduced in the data to avoid the effects of this deviations. Moreover we have 
added the simplest type of correction term which arises from the replacement 



we use ( |1.3|) . 7 does not appear as a subleading correction to fi c and one can 
use canonical Monte Carlo simulations (the so-called link-flip algorithm ||) which 




(2.1) 
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in (JO) and fitted to the form 

\n(n NT ) = A + (7 - 2) ln(2J(l - jL)) + - (2.2) 

for B > B . A and C are some fit parameters. Comparison of the results with 
and without this correction term can be seen in fig. 2 where we plot the value of 7 
extracted with different cut-off's B . We see that including the correction improves 
the results considerable. 

Let us assume that the values jb extracted from (|2.2| ) appproach exponentially 
a limiting value for large B : 

7s =7- Cl e- C2B °. (2.3) 

The result of such a fit is shown in fig. 2. It is clear from fig. 2 that the assumption of 
an exponential approach of 7s to 7 is not essential for the extraction of 7. We have 
introduced it at this point in order to treat all measurements consistently. For the 
matter fields coupled to gravity the finite size effects will be larger and extrapolation 
to large Bq more important. 

The 7 extracted in this way for different lattice sizes is: 



Nt 7 

1000 -0.496 ± 0.005 

2000 -0.501 ±0.004 

3000 -0.504 ±0.004 



which is in good agreement with the expected value of 7 = —0.5. It shows that this 
kind of simulations are indeed well suited to measure 7 and it is thus natural to try 
apply them to the case of matter couple to 2d gravity. 



2.2 The Ising model 

The next non-trivial test of the method is to study the Ising model coupled to 2d 
gravity. It has been solved analytically || and was found to have a 3r<i-order phase 
transition. The coupling to gravity is in a sense weak as it only changes the string 
susceptibility at the critical point (from 7 = —1/2 to 7 = —1/3). For this reason it 
has until now been considered very difficult to measure 7 directly, since it required a 
fine-tuning of both the bare cosmological constant fi and the spin coupling constant 
(3. On the other hand it has been verified that it is indeed possible to extract the 
other known critical exponents [^] since for these exponents it is possible to use the 
canonical ensemble in the simulations. 
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The Ising spins are placed in the center of the triangles and they interact with 
the spins on neighbouring triangles. This corresponds to placing them on vertices 
in the dual graph. In that case the critical point has been found explicitly and is 
(3 C = 0.7733... The (canonical) partition function of the model is 



Zn t (P)= E I' ^ (2-4) 

T~N T to} 

where the summation is over all triangulations with Nt triangles. 



In the simulations we used a Swendsen-Wang cluster algorithm |12|] to update 
the Ising spins and lattices sizes Nt = 1000 and 2000. We made runs for several 
values of the coupling in the interval 0.6 < (3 < 0.95 and then fitted the distributions 
to eq. ( p.2|) . In this way we could extract values 7b (/3) and by assuming a relation 
like (|2j): 

JbM = >y<J3) - Cl ((3)e- cmBo (2.5) 

we have extracted the values for 7(/3) shown in fig. 3. Examples of 7_b (/3) and 
the exponential fit ( |2.5| ) for different values of (3 are shown in fig. 4. We observe a 
marked increase in the dependence on B when j3 approaches (3 C . 

We get, as expected, the pure gravity value of —0.5 for couplings far below and 
above f3 c . In the vicinity of the phase transition we see on the other hand a clear 
peak and the peak values agree well with the exact value 7 = —1/3. We conclude 
that the method for extracting 7 works well in this case too, although it should be 
clear that the amount of numerical work needed is much larger in this case than in 
the case of pure gravity. 



2.3 Gaussian fields 

The gaussian fields x^, fi = 1, . . . , D are placed on the sites i of the triangulation 
T. They can be viewed as representing an immersion i — > xf of our abstract tri- 
angulation T into R D , i.e. a model for non-critical strings and they also represent 
a coupling of matter with central charge c = D to gravity. The multiple gaus- 
sian models do not interact directly with each other but only through their mutual 
interaction with the geometry. The (canonical) partition function is given by 

z *t= e / n '/'.'•,< - • , ' ;:: * (2.6) 

where the summation is over all triangulations T with Nt triangles. One site is kept 
fixed in order to eliminate the translation mode. No coupling constant appears in 
the action as it can be absorbed in a redefinition of the gaussian fields. 
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Again we have performed simulations with up to 10 7 sweeps for lattice sizes 
ranging from 1000 to 4000 triangles. We have used from one to five gaussian fields 
and a standard Metropolis algorithm to update them. In fig. 5 we show how 
the distributions of baby universes change with increased c (normalized with the 
distribution for pure gravity). Fitting these distributions to the functional form 
Q2.2|) and extracting 7 as above yields the results shown in fig. 6. The results are 
compatible with earlier estimates ||. 

It is seen that 7 is too small for c = 1 where it is known that 7 = 0. But in the 
case c = 1 we know that the asymptotic form is not correct. It should be 



multiplied with logarithmic corrections. If we include these we get for c = 1 that 
(O) is replaced by 



n NT (B) ~ [(N T - B)B} 1 ~ 2 [ln(iV T - B) hiBf . (2.7) 



In this formula we have left 7 and a as variables. Model calculations give a = —2, 
but it is not known whether this power is universal and the model has not been 
solved analytically in the case of one gaussian field. 

If we fit to ( |2.7p in the way described above (including also the 1/B correction) 
we extract for c = 1 the following values of 7 and a for different lattice sizes: 

Nt 7 a 



1000 -0.22 ±0.05 -0.5 ±0.4 
2000 -0.14 ±0.07 -1.0 ±0.4 
4000 -0.09 ±0.08 -1.2 ±0.4 



Both 7 and a moves towards the expected values and -2 as a function of Nt, but 
the finite size effects are clearly larger here than for pure gravity. 

In fig. 6 the results of a fit to Q2.7D for c > 1 is included. It is seen that 7 
extracted in this way exceeds the theoretical upper bound 7 = 1/2(0). In addition 
the power a decreases from -1.2 for c = 1 to -5 for c = 5. We conclude that either 
logaritmic corrections are not the right ones to include for c > 1 or finite size effects 
are so large that they make the fits unreliable. 

What is clear from the analysis is that 7 increases with c for c in the range — 5. 
According to ( |1.3|) this means that the number of baby universes of a given size will 
increase, i.e. the fractal structure will be more pronounced with increasing c. We 
have illustrated this in fig. 7, which shows two "typical" surfaces corresponding to 
c = and c = 5. It should be emphasized that the pictures are only intended to 
visualize the internal structure, i.e. the connectivity of the surfaces[| 

lr The surfaces are constructed in the following way: Given the connectivity matrix of the trian- 
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3 Discussion 



We have verified that the technique of extracting the entropy exponent 7 directly 
from the distribution of baby universes is superior to the methods used until now 
from a practical point of view. A single (although long) Monte Carlo run for a fixed 
value of Nt is sufficient for extracting 7 and we get the correct results for c < 1. 
On the other hand the situation in the case c > 1 has not really improved much 
compared to the earlier measurements ||. We get in fact similar results, and this 
shows that the method also works in the case c > 1 and the ambiguity in extracting 
7 for c > 1 is that we do not know the correct functional form to be used in the 
fits. It is clear that it would be most interesting if we could reverse the procedure 
and use the data to obtain knowledge about the corrections to ( |L1|) for c > 1. Our 
data are not yet good enough to do this in a convincing way, but the problem is 
clearly not due to the baby universe technique introduced in this paper, but due to 
the inefficiency of the flip algorithm used to update the triangulations. 
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Figure captions 

Fig.l The distribution of baby universes in the case of pure gravity. 

Fig. 2 Fitted values of 7s for different cutoff's B for pure gravity. Values are shown 
for fits with and without the correction term included. The curve shows a fit 
using ( |2.3| ) resulting in 7 = —0.496 ± 0.005 (errors are 95% confidence limits 
for a x 2 -test). 

Fig. 3 Fitted values of 7 vs the coupling (3 in the case of one Ising model coupled to 
gravity. Results are shown for two lattice sizes, N T = 1000 and 2000. 

Fig. 4 Fitted values of Jb {P) f° r various (3 as a function of the cut-off B . The curves 
represent fits to (|2.5| ). 

Fig. 5 The distributions of baby universes for up to five Gaussian fields coupled to 
2d gravity. The values are normalized with the distribution for pure gravity. 

Fig. 6 Fitted values of 7 vs central charge in the case of multble Gaussian fields. 
Results are shown for fits without (circles) and with (squares) a logarithmic 
correction term included. 

Fig. 7 3d illustration of the fractal structure of the surfaces for c = (fig. 7a) and 
c = 5 (fig. 7b). N T = 200 is used. 
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